Momentum distribution of liquid helium 
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We have obtained the one-body density matrix and the momentum distribution n(p) of hquid *He 
at T = 0°K from Diffusion Monte Carlo (DMC) simulations, using trial functions optimized via 
the Euler Monte Carlo (EMC) method. We find a condensate fraction smaller than in previous 
calculations. Though we do not explicitly include long-range correlations in our calculations, we get 
a momentum distribution at long wavelength which is compatible with the presence of long-range 
correlations in the exact wave function. We have also studied "^He, using fixed-node DMC, with 
nodes and trial functions provided by the EMC. In particular, we analyze the momentum distribution 
n(p) with respect to the discontinuity Z as well as the singular behavior, at the Fermi surface. We 
also show that an approximate factorization of the one-body density matrix p(r) ~ po{r)pB{r) holds, 
with po{r) and pB(r) respectively the density matrix of the ideal Fermi gas and the density matrix 
of a Boss ^He. 

PACS numbers: 67., 67.40.-w, 67.55.-s 



I. INTRODUCTION 



:< ajjOp > is a fun- 



The momentum distribution n(p) 
^Ciiamental quantity for the study of both the static and 
'J^he dynamical properties of quantum liquids, as it gives 
Q(iirect information on the high rnomentum components 
|t)£ the ground state wave functioru. Experimentally, ac- 
'^ess to n{p) is provided by deep inelastic neutron scat- 
C3;ering at large momentum transfer TiQ. The extraction 
Q)f n(p) from the measured scattered intensity, however, 
Hs] affected by the the limitations imposed by the experi- 
K*lnental resolution and the final-state interactions. Thus, 
kJfhe most accurate information on the momentum distri- 
*Vj>ution of "^He available to date is likely to be the one ob- 
Cd^ained through accurate, microscopic calculations, such 
as those presented in this study. 

At T = 0°K the momentum distribution n{p) of an 
ideal Bose gas is given by a delta function (5(p), corre- 
sponding to all particle being in the condensate. On the 
contrary, the condensate fraction of "^Hc at the equilib- 
rium density is less than 10%, implying that the effects of 
the strongly repulsive core of the inter-atomic interaction 
is non perturbative. Similarly for "^He the discontinuity 
of n{jp) at the Fermi momentum fipp, which gives the 
strength of the quasiparticle pole, is ^ 0.2 at equilibrium 
instead of 1, as in the ideal Fermi gas. In fact, it has 
been always difficult to carry out ah initio calculations 
of the momentum distribution, for these systems, within 
the field-theoretical approacho. 

Modern, realistic, quantitative calculations only 
started with the development of Variational Monte 
Carlo (VMC) methodifLpjid the Hyper-Netted-Chain 
(HNC,FHNC) equationsaEl for Jastrow models of both 
Bose and Fermi liquids. FHNC and HNC equations 
for the momentum distribution have been numerically 



solved iEui/ariational calculaijaas of nuclear mattem, liq- 
uid ''HeQl^, and liquid ^HeItEI, for which triplet and 
backflow correlations were also taken into account. To 
improve upon the above variational estimates of n[p), 
non conventional pertiwbative techniques, based on cor- 
related basis functionsO|JCBF), have been i^veloped, 
and applied to liquid '^HelZl and nuclear mattem-El. 

HNC and FHNC theories have the merit to allow for 
fine details of interparticle correlations, such as long- 
range behavior, spin dependence, anisotropics in inho- 
mogeneous systems. However, not all cluster diagrams 
resulting from the theory and involved in the HNC for- 
malism, can be summed in closed form and in FHNC the 
procedure to estimate the elementary diagrams with ex- 
change bonds is not completels^ under control. |-|Therefore, 
approximations like. .sealingiS, interpolationa or trun- 



cated summations 



must be invoked, and, in liq- 



uid helium calculations, these result in a non negligi- 
ble loss oi-actjuracy with respect to a fully variational 
treatmentll3'LD. 

The momentum distribution of liquid ^He and '^He has 
also been calculated, at zero temperature, by ijswg the 
Green Function Monte Carlo (GFMC) methodlatt At 
finite temperature calculations have been jptjcformcd for 
"He by Path Integral Monte Carlo (PIMC)HEI and, more 
recently, by VMC with trial functions of the Shadow 
typec3. The available theoretical estimates of the mo- 
mentum distribution of liquid ^He provided by varia- 
tional and GFMC methods are in reasonably good agree- 
ment among themselves, except for low momenta and for 
the condensate fraction. For liquid ^He the situation is 
less satisfactory. 

The GBMCUand the Diffusion Monte Carlo (DMC) 
methodg2i'c2[ll3, afford to date the most precise tools 
to perform ground-state calculations for many-particle 
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systems. For Bosons they provide estimates of the en- 
ergy that are virtually exact, within the statistical accu- 
racy. In fact the same is true for the averages of oper- 
ators that are (R-)diagonal (not the case for n{p)\), for 
which algoi jit | hn is exist, like the so-called forward walk- 
ing methodcaa, that yield pure estimates. On the other 
hand, the estimates of observables that do not commute 
with the Hamiltonian are usually obtained from mixed 
averages — through an extrapolation procedure whose ac- 
curacy depend on the qualitynof the trial function used 
for the importance samplingH. The extrapolation intro- 
duces a bias in the estimates which is second order in the 
difference between VE* and ground-state wave function $0. 
For Fermions there is an additional source of error related 
to the existence of the so-called sign problem. To date, 
to get a numerically stable algorithm it is customary to 
approximate the unknown nodes of the sought ground 
state $0 with jLhose of the trial function ( fixed-node 
approximationE3). This imposes a bias on any average, 
which for non-diagonal observables cumulates with the 
one arising from the extrapolation procedure. Therefore, 
to minimize systematic errors, especially in the evalua- 
tion of properties such as the momentum distribution, it 
is necessary to achieve maximum accuracy in the opti- 
mization of the trial function. 

Recently, a new optimization procedure based on 
Monte Carlo calculations and denoted pas-JCuler Monte 
Carlo (EMC) method has been proposedE§llZl. This EMC 
method iiaa-been successfully applied to both liquid "'He 
and ^HelijO. The EMC wave functions have pair and 
triplet correlations fully optimized and provide the low- 
est available energy upperbounds. Moreoveij, i thftirj use 
in Diffusion Monte Carlo (DMC) calculationsHiNllZl has 
led to results of unprecedented accuracy for the energy, 
pair function and static structure function. 

In this paper we present results for the one-body den- 
sity matrix p(r) and the momentum distribution n{p) of 
liquid ''He and ^He, at various densities, obtained with 
DMC calculations based on EMC wave functions. The 
plan of the paper is as follows. In the next Section we 
summarize the computational details involved in the cal- 
culation of n{p). We then present the results for liquid 
"^H in Sec. Ill and those for liquid ^He in Sec. IV. We 
finally offer a summary and conclusions in Sec. V. 

II. COMPUTATIONAL DETAILS 

In DMC simulationgiil@ the ground state wave func- 
tions $0 is sampled through a random walk in configu- 
ration space, guided by a trial function 5". In practice 
one samples the mixed probability / — <I>oVE'. We chose 
^ = SF for ''He and * = D^DiF for ^He, with the cor- 
relation part F symmetric in the particle coordinates, S 
a symmetrized product of one-particle orbitals and D| 
and Di Slater determinants of one-particle orbitals for 
particles of up and down spin. In the homogeneous liq- 
uid S = const and D„ is built from plane waves (PW), or 



from plane waves with short (SBF) or long (LBF) ranged 
back-fiow correctjans. As a full account has already been 
given elsewheretj'EZl of both the EMC method, which we 
employ to construct and optimize the trial wave func- 
tion \E', and of the use of EMC wave functions in DMC 
simulations, here we shall restrict to essential details. 

All the calculations (variational and diffusion) pre- 
sented in this paper, unless explicitly noted, have been 
performed using EMC wave functions with fully opti- 
mized pair and triplet pseudopotentialsllZl (OJOT)and 
modeling pHelium with the HFDHE2 pair potential of 
Aziz et a/.Ej. A cubic simulation box and periodic bound- 
ary conditions werp-used. For ^He bapkfktw correlations 
have been includedllZl, in the usual wayEZrEj, by replacing 
the plane waves exp (ik^ • Tj) in the Slater determinants 
with exp {iki ■ Xj), where = rj+J2k=ij vifjk^j -ffc)- 
The function ri{r) is taken either short-rangecO (SR) 

ris{r) = \b exp (-(r - rsf/^l) ((2r - L)/L)\ (1) 

or long-ranged0 (LR) 

77i(r) = \b exp (-(r - rsf/^l) + A^j A^ (2) 

with As, rs, and A'^ variational parameters. The 
long-ranged backflow function is smoothly cutoff at the 
the boundary of the simulation box of size L by replacing 
the expression given in Eq. (^) with rj'j^{r) = rj^ir) -\- 
rjL{L — r)~2r]L{L/2). In practice, we first simultaneously 
optimize the pair and triplet pseudopotentials with the 
backflow parameters as specified in Refs. ^ and ^ and 
then we optimize ri{r) at fixed pseudopotentials. 

The "^He results presented below were obtained with 
short-ranged backfiow, unless otherwise specified. Also 
the DMC simulations were performed within the fixed- 
node approximation, whereby the nodes of the ground- 
state $0 are assumed to coincide with those of the EMC 
trial function. 

In a uniform liquid in a state described by the wave 
function ^ the one-body density matrix can be defined 
as 

p{r) = 

V J rfr2 ■ ■ ■drAf^*(ri,r2, ■ ■ ■ ,rAf)^(ri + r, r2, • ■ ■ , rjy) 
/ dri • • •drAr|^'(ri,r2, • • • ,rAr)|2 

(3) 

so that, having imposed periodic boundary conditions to 
the N particles in the volume V, 

P(0) = 1, (4) 

and the independence on ri follows from translational in- 
variance. For Fermions, the integration over the variable 
Ti is understood to imply also the trace over spin projec- 
tion. If we denote with R and R' respectively the config- 
urations (ri, r2, • ■ • , rjv) and (ri -I- r, r2, • • ■ , tn), and we 
exploit the translational invariance, we can rewrite Eq. 
i) as 
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p{r) 



J dm*{R)-i'{R') 



dRP{R) 



with 



P{R) = 



JdRl^iRW 



(5) 



(6) 



the probabiUty induced by the wave function 'J. Us- 
ing Eq. the variational density matrix (i.e., the one 
defined in term of the trial wave function 5') may be 
conveniently calculated by Monte Carlo as 
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(7) 



with the configurations Ri drawn from the probability 
P{R)- The DMC extrapolated estimate p{r) is given by 



p{r) = 2/9„iix(r) - pvMc{r), 



(8) 



where the mixed estimate Pinix('') is calculated from an 
expression identical with that of Eq. (|^), with the con- 
figurations Ri drawn however from the mixed probabil- 
ity PmUR) = fiR)/ JdRfiR), /(i?) = $o(i?)1'(/^). In 
practice, the auxiliary configurations R' appearing in Eq. 
(|^) are generated from a given R by moving a particle ei- 
ther of fixed increments along a random direction (FM) 
or to points randomly distributed in the simulation box 
(RM). It turns out that the two methods give more ac- 
curate results at small and at large r, respectively. 

The momentum distribution is defined as n(j)) — 
(ajjOp), where an average on spin projections is also im- 
plied for unpolarized "^He. It is simply related to the 
one-body density matrix, by a Fourier transform: 



n{p) = ^ J dre*P >(r) 



p) + 1 dre^P-^ {p{r) - no)) , 



(9) 



with p — N/V the density of the system, f the degen- 
eracy factor, which is 1 for ''He and the fully polarized 
•^He, and 2 for normal ^He, and uq the large r limit of 
the density matrix, uq — p{oo). In fact p(r) vanishes for 
large values of r in ^He, whereas for ''He it saturates to 
riQ ^ 0, due to macroscopic occupation of the state with 
zero momentum, uq is the condensate fraction, i.e., the 
fraction of ^He atoms occupying the state with p = 0. 

Evidently, the normalization of the density matrix 
given in Eq. (^) implies, for the momentum distribution, 
the normalization sum rule 



(27r)3p 



dp nip) =p(0) = l. 



(10) 



The momentum distribution has been calculated in two 
different manners. Having sampled the density matrix 
as function of r one can just take its Fourier transform 



according to Eq. (|9|). Alternatively, restricting to the 
RM method, one can also directly accumulate 
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^*(i?0 



(11) 



At a first sight it seems that the r integration is missing 
in Eq. (pTl) above. However, a little reflection shows that 



accumulating the estimator of Eq. (11) correctly imple- 
ments the integration (average) over r, which appears in 
the definition (^) of the momentum distribution, since 
r is chosen at random for each configuration Ri, with a 
uniform distribution in the simulation box. This second 
calculation of n(p) is implemented for p's that are re- 
ciprocal lattice vectors of the simulation cell. Eq. ( pi] ) 
is employed to accumulate variational and mixed estima- 
tors, from which the extrapolated estimator is then ob- 
tained, as explained above. From Eqs. it follows that 
the kinetic energy per particle T can be related to the 
curvature of the density matrix at the origin, according 
to 



T = - 



(12) 



r=0 



which also implies the kinetic energy sum rule (see Eq. 



2M (27r)3p 



dpp^n{p). 



(13) 



III. LIQUID *HE 

We have carried out DMC calculations of the one-body 
density matrix and of the momentum distribution of "^He 
at four densities, using EMC trial functions and 64 atoms 
in the simulation box. Selected runs with up to 232 par- 
ticles have been performed to check for finite size effects. 
In Fig. 1^ we show our results for the one-body den- 
sity matrix p(r) at the equilibrium density. The kinetic 
energy sum rule ( p^ ) is manifestly satisfied and the sat- 
uration to a finite no at large r is evident, in spite of the 
fact that with 64 particles only distances up to about 7 A 
are accessible. We also give in the inset a comparison 
between estimates obtained with the FM and RM meth- 
ods. The greater accuracy of the latter method at large 
distances is apparent. In Fig. ^ we report, also at the 
equilibrium density, extrapolated estimates of the mo- 
mentum distribution n(p), obtained using Eq. (pT|). Due 
to the finite size of the system only wave vectors larger 
than ~ 0.4A~^ are accessible. We note that a shoulder 
is discernible in pn{p) at p > 2A. 

In order to extract the condensate fraction riQ from 
our DMC results, as well as to facilitate applications, we 
have fitted our DMC extrapolated estimates of the den- 
sity matrix and momentum distribution, obtained with 
the RM method, to the following analytic formula: 
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FIG. 1. One body-density matrix of He at the equilibrium 
density p(A"^) = 0.02186. The full curve is the fit of Eq. Q 
to our DMC results (extrapolated estimates) , with the dashed 
curve showing the parabola 1 — {MT/3Tt^)r^ that satisfies 
the kinetic energy sum rule. The solid circles and rectangles 
respectively give our DMC results obtained with the FM and 
RM methods. 
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The first 2 terms in Eq. (|4|) account for the existence 
of the condensate, while the third suitably models the 
shoulder in pn{p). The fourth and last term accounts for 
the gross main structure of the momentum distribution. 

We have simultaneously fitted nij)) and p(r), which 
must be obtained numerically by 

Fourier transforming the function of Eg . ( [1^ ) , imposing 
as well the normalization condition ( |lO| ) and the kinetic 
energy sum rule (p^). Moreover, we set pi — Mcjlh, to 
satisfy the long-wavelength behavior 

TABLE I. Parameters of the fit (|l|) to the DMC momen- 
tum distribution and density matrix of ''He at T = 0°K, at 
various densities, p is in A""' and pi-ps arc in X^^ . 



p 


0.01964 


0.02186 


0.02401 


0.02622 


no 


0.11163 


0.071673 


0.046227 


0.027079 


a 


1.6941 


1.7634 


2.2342 


1.9976 


ni 


0.055274 


0.038282 


0.019157 


0.021497 


712 


0.40285 


0.39893 


0.36695 


0.38483 


Pi 


0.49957 


0.73394 


0.95934 


1.1964 


P2 


0.38971 


0.38525 


0.40435 


0.40588 


P3 


1.4900 


1.6744 


2.2153 


2.1832 


P4 


0.29713 


0.20538 


0.28494 


0.062765 


Ps 


0.85128 


1.0016 


1.0132 


1.2549 




FIG. 2. Momentum distribution of *He at T = 0°K. Long 
dashes, full curve and short dashes are fits to the DMC results 
at p(l"^) = 0.01964, 0.02186, 0.02622, whereas the rectangles 
give the DMC extrapolated estimates at the equilibrium den- 
sity p{A~^) = 0.02186. The inset shows the long-wavelength 
behavior pi n(p), with the dotted curve reporting the results 
of GFMCB. 



lim pn{p) — 

p—^Q 



HqMc 

2h ' 



(15) 



induced by long-ranse,correlations, as first discussed by 
Gavoret and NozieresEd. Above M is the mass of the ^He 
atom and-c the sound velocity, which we estimate from 
the DMCO equation of state (EOS). Thus we fit at each 
density 6 independent parameters to more than 100 MC 
points, obtaining a reduced between 0.98 and 1.20. 
The resulting fit parameters are recorded in Table ||. 

In principle one could take c too as unknown and get 
an independent estimate of the sound velocity. We have 
tried this alternative, at the equilibrium density, obtain- 
ing an estimate of c that is about 10% lower than ex- 
perimental and DMC-EOS estimates, and has however a 
very large uncertainty ( ~ 50%), reflecting the absence 
of DMC points for p < 0.4A~-'^), where the singular term 
in n{p) is important. We may conclude that our data a*e 
compatible with the presence of singular term in n(p)l22l, 
in spite of the lack of long-Mnge terms in the pseudopo- 
tentials that we have usedllll. In Fig. ^ we show the 
density matrix at the four densities that we have stud- 
ied, as given by the fit of Eq. (|l^) with the parameters 
of Table |. 

DMC momentum distribution and density matrix are 
compared with the fit of Eq. (14) in Figs. |l|, ^, and ^, at 
the equilibrium density. The fit appears to be very good. 
A similar conclusion holds at the other densities that we 
have studied. 

The condensate fraction hq is mostly constrained by 
the large r behavior of p(r), which results into the term 
proportional to <^(p) in the momentum distribution. As 
we have already mentioned, the singular behavior of nij)) 
at small p, implied by Eq. (p^), is much less effective in 



4 




0.12 



0.08 



0.04 



0.00 



FIG. 3. One-body density matrix p{r) of ''He. The curves 
from the topmost to the lowest, give the analytical fit of Eq. 
jl^ ) to our DMC results (extrapolated estimates), respec- 
tively at p(A"^) = 0.01964, 0.02186, 0.02401 and 0.02622. 
The inset shows p{r) in the tail region with an enlarged scale. 
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FIG. 4. Dependence on the quality of the wave func- 
tion and on the method (VMC versus DMC) of the one 
body-density matrix of ''He, at the equilibrium density 
p(l"^) = 0.02186. The full curve is the fit to our DMC re- 
sults (extrapolated estimates) for 64 particles, using our best 
trial function (OJOT), while the dashed curve gives the fit 
to the VMC results obtained with this trial function. Empty 
rectangles and solid squares give respectively VMC and DMC 
results obtained from a simple OJ trial function, which em- 
bodies only pair pseudopotentials. Finally, the triangles are 
the DMC results obtained using the OJOT trial function and 
232 particlasn and the circles report the finding of PIMC at 
T= 1.18°kI13. 

determining no, because of the absence of DMC estimates 
for p < 0AA~^. From the inset of Fig (|^) one might 



conclude that the range r ^ 7 A accessible with 64 atoms 
is not big enough too unambiguously assess the value of 
no- However, simulations with 232 atoms yield, in the 
extra range 7 A ^ r ^ 11, DMC estimates that are in 
perfect agreement with the fit to the 64 particles results, 
as it is clear from Fig. ^. We also illustrate in this figure 
the dependence of the large r limit of density matrix, i.e., 
TiQ, on the quality of the wave function. In particular it 
is apparent than improving the MC description cither 
changing from VMC to DMC, for given trial function, or 
changing to a better trial function in DMC, results into 
a decrease of no, in the case considered. 
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FIG. 5. Condensate fraction of liquid 4He, as a func- 
tion of the deasity. DMC (solid cij-cles); EMC (open cir- 
cles); GFMCta (diamonds); JJNCQ (crosses); experimen- 
tal estimates, at T = 0.75°KE3 (solid triangles); PIMC at 
T = 1.18°Ky (empty triangle). 

In Fig. I and Table || we compare our predictions for 
no with those from other theoretical treatments as well 
as with experimental results at low temperature. Consis- 
tently with the observation made above, our use of very 
accurate trial functionstll yields DMC predictions for the 
condensate fractioa which areJower than previously ob- 
tained by GFMCta and HNCQ. On the theoretical side 
the only prediction that agrees with our own, though it 
has a much larger statistical error, is the PIMC one at 
T — 1.18°KI13. We should remind the reader that in 
fact PIMC has no trial function bias. We find instead a 
sizeable discrepancy with the experimental estimates of 
Snow et a/L^I, who determine no by fitting a model n{p) to 
the measured Compton profile J{Y) — a procedure how- 
ever which appears to be model dependent. Differeat 
choices for n{p) produce equivalently good fits of J(y)E3, 
though embodying very different condensate fractions, all 
the way from no = to no = 10%. In fact , our DMC 
momentum distribution yields prediction for the Comp- 
ton profile which agrees well with the experiments, as we 
show below. 
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FIG. 6. Compton profile J{Y) of *He at tlie equilibrium 
density p{A~^) = 0.02186 compared with exp|e*imental (solid 
circles) data at Q = 23 A'^ and T = 0.35°-ft£l. The dashed 
(solid) curve is obtained from our fit to the DMC n{p) allowing 
(not allowing) for both the experimental broadening and the 
final state interactions of Ref. |^ In calculating the dashed 
curve a shift of —0.1 A in Y has been also used as in 

The inelastic neutron scattering cross-section at high 
momentum transfer HQ can be approximated by its Im- 
pulse Approximation (lA) expression, which is propor- 
tional to the Compton profiled 

J{Y)^^J^^^dppn{p). (16) 

The scattering, in lA, does not depend on the energy 
Lo and the momentum transfer TiQ separately, but only 
through the scaling variable Y ^ given by 

Y ^ {M/hQ){uj ~ujr), (17) 

where lo^ — fi'^Q^ /2M is the recoil energy of the scat- 
tering atom. The dynamical response function S{Q,uj), 
in lA, is given by J(Y) times the factor M/{hQ). Final 
state effects (FSE) of the medium on the scattered atom 
as well 

TAB*,E II. Conjdensate fraction no in "He. DMG, 
GFMCta, and HNCfl predictions are at T = 0. The PIMCta 
result is at r = 1.18°K and the density is in 1-=^. The figure 
in parenthesis is the uncertainty on the last figure, whenever 
available. 



p 


0.01964 


0.02186 


0.02401 


0.02622 


DMC 


0.112(1) 


0.0717(5) 


0.0462(6) 


0.02.71(6) 


PIMC 




0.069(10) 






GFMC 




0.092(1) 


0.052(1) 


0.037(2) 


HNC 




0.092 


0.065 


0.043 



as experimental resolution (ER) broaden up the Comp- 
ton profile, particularly its delta peak at p = 0, which 
is due to the Bose condensation. In Fig. || Compton 
profiles, calculated with our DMC momentum distdbu- 
tions are compared with observed scattering dataO at 
T = 0.35°K, convettfid to J{Y). Once ER and FSE are 
taken into accounto, good agreement with the experi- 
ment is obtained . 

IV. LIQUID ^'HE 

For normal '^He, DMC simulations using EMC trial 
functions with backflow and the fixed-node approxima- 
tion have been performed at five densities, with 54 atoms 
in the simulation box. We have investigated the depen- 
dence of the momentum distribution on the size of the 
system and on the range of the backflow, respectively 
with runs for 114 atoms and with runs using trial func- 
tions embodying long-range backflow. 
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FIG. 7. One body-density matrix of ^He at the equilibrium 
density p(A"^) = 0.01635. The full curve is the fit of Eq. (^ 
to our DMC results (extrapolated estimates) , with the dashed 
curve showing the parabola 1 — {MT/3h'^)r^ that satisfies 
the kinetic energy sum rule. The solid circles and rectangles 
respectively give our DMC results obtained with the FM and 
RM methods. 

In Fig. 1^ we give our DMC estimates for the density 
matrix at the equilibrium density. It is clear that the 
kinetic energy sum rule ( |l^ ) is satisfied. The size of the 
system allows for the determination of p{r) through its 
first zero and up to the first minimum. A comparison 
between results obtained with FM and RM methods is 
also given, in the inset. The smaller error on the FM 
results, compared with that on the RM estimates, is due 
to the much longer runs used to accumulate the FM p{r) 
in this case. 
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FIG. 8. Momentum distribution of ^He at T = 0°K 
and p(l"^) = 0.01413, 0.01635 and 0.01946: dots are DMC 
(fixed-node) extrapolated estimates; full curves are fits to the 
DMC results. Larger densities correspond to lower curves 
at p = 0. The inset shows the behavior around the Fermi 
momentum at the equilibrium density p(A~'^) — 0.01635A^, 
using PW (• ■ •), SBF (— ) , and LBF ( ) trial functions. 

In Fig. H we report, at three different densities, the 
momentum distribution n{p) obtained using Eq. (|ll|). 
Clearly, the discontinuity Z at the Fermi wave vector pp 
is substantially reduced, with respect to its value Z = 1 
in the noninteracting system, and moreover it systemati- 
cally shrinks, as the system gets denser and the effects of 
the interaction become more important. Z is also slightly 
reduced 




FIG. 9. Momentum distribution of He at equilibrium den- 



sity p{A 



0.01635 with the VMC method, using OJ trial 



functions, with 54 (solid circles) and 114 (rectangles) parti- 
cles. 



when the nodes of the trial function are improved from 
PW to SRB and then to LBF. Size effects on the momen- 
tum distribution appear to be negligible, as it is clear 
from the comparison between variational results for 54 
and 114 particles given in Fig. ^, at the equilibrium den- 
sity. 

In a normal Fermi liquid such as '^He the momentum 
distrihjitkin, in addition to the discontinuity, has infinite 
slopea^j'EEl at pf- To leading order in p — pp 



PF 



In 



Pp 



P-PF 



Pp 



(18) 



with the coefficient A related to imaginary part of the 
self-energy i;(p, E). 

To extract Z from the calculated momentum distribu- 
tion and to check that our results are consistent with the 
presence of the singular term of Eq. ( |l8| ) we have fitted 
our data for p(r) and n{p) to the real space form 



p{r) = Zpo{x) + a 



Pa{x) xpo{x) — sin{x) 



02- 



-ai- 



+ (03 + a4X + a^x^ + aQX^)e 



with X = ppr, and 



Po(^) — -^(sinx — xcosx), 
x-^ 



(19) 



(20) 



the density matrix of the ideal Fermi Gas. The first two 
terms in Eq. (^9|) account respectively for the discon- 
tinuity and the infinite slope of n{p) at pp. The third 
term allows for a finite discontinuities in the first and sec- 
ond derivatives of n{p) at pp. Finally, the fourth term is 
needed to eliminate the divergence that the term po{x)/x 
produces at the origin. We impose the normalization con- 
dition (|l^) and the kinetic energy sum rule ([l^), as well 
as the vanishing of the first and third derivative of p(r) 
at r = 0. Thus we fit at each density 5 independent pa- 
rameters to more than 100 MC points, with a reduced 
between 0.87 and 1.15. The resulting fit parameters are 



recorded in Table III 



We compare fit and DMC estimates for p{r) and for 
n(p), respectively in Fig. ^ and Fig. ||. Evidently, Eg. 
( |l9| ) is fully consistent with our DMC results. In Fig. |o| 
we show the one-body density matrix at three of the den- 
sities that we have studied, as given by the fit of Eq. (|l|). 
As for the non-interacting case, p(r) becomes steeper 
with increasing the density and its first zero moves to- 
ward the origin. 

A comparison of our results for the momentum distri- 
bution of '^He with those from some other calculations 



is given in Fig. 11 and in Table IV. HNC resultgl3 are 



in close agreement with our variational n(p), for p > pp 
as well as with the estimate of Z, while small differences 
are present at small momenta, which are however of little 
relevance in the density of states oc n{p)p^. The evidait 
discrepancies between our results and those of GFMCEZI 
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FIG. 10. One-body density matrices p of ^He. The curves 
from the less to the most steep, give the analytical fit of Eq. 
( |l9| ) to our fixed- node DMC results (extrapolated estimates), 
respectively at p(A"^) = 0.01198, 0.01635 and 0.01946.The 
inset shows p(r) in the tail region with an enlarged scale. 

around the Fermi surface are pcobably due to the poorer 
trial function used in GFMCEZI, as well as to statistical 
errors. 

As we have already mentioned, the functional form of 
Eg. ( ]l9| ) implies for n{p) a singular term sipp of the form 
dig), with A — ai/n. For instance, at the equilibrium 
density p{)r^) = 0.01635 we find A = 0.06(2), being 
however unable at present to assess the size dependence 
of such estimate. An independjent estimate of A is given 
by a perturbation calculationE3 for a dilute, hard-sphere, 
Fermi gas. To order {ppR)'^, 

Ao = 2^^{ppRf, (21) 

with R the radius of the Fermi particle. Evidently '^He 
at equilibrium is no at all dilute. Nevertheless, taking 
R = 0-/2 ~ 1.3l, one gets - 0.20, which is of the 
same order of magnitude as our DMC estimate. In CBF 



TABLE III. Parameters of the fit (|l^ to the fixed-node 
DMC momentum distribution and density matrix of '^He at 
T = 0°K, at various densities, p is in A"''. 



p 


0.01198 


0.01413 


0.01635 


0.01797 


0.01946 


z 


0.45977 


0.30598 


0.23616 


0.14328 


0.13566 


ai 


0.14248 


0.28612 


0.19376 


0.41699 


0.24613 


a,2 


-0.11973 


-0.21912 


-0.12381 


-0.28125 


-0.23982 


13 


-0.73200 


-0.80732 


-0.64305 


-2.1780 


-0.96272 


a4 


3.2095 


1.9033 


3.2723 


4.1029 


3.9981 


as 


-1.5905 


0.64801 


-0.098765 


-2.4600 


-0.64725 


aa 


7.6229 


1.5582 


4.6639 


6.8026 


6.2113 


bi 


9.0415 


5.3494 


7.3460 


7.3676 


7.5530 


b2 


3.5907 


2.7487 


3.0716 


3.3316 


3.1652 



0.6 




0.5 1 1.5 2 



p a-') 

FIG. 11. Momentum distribution of ^He at equilibrium 
density p(l"^) = 0.01635: comparison of the present 
DMC Uy) and VMC (- r-,-) fits with the predictions of 
GFMCy(circles) and HNcta (solid circles). 

theory one obtainsi, instead, 

_2Wo de^ipp) 

AcBF — Pf ; , y^^) 

■nap 

where e"{p) is the variational energy and Wg is an in- 
verse energy parameter characterizing the imaginary part 
W{p, E) of the self-energy I](p, E), close to the Fermi en- 
ergy, 

W{p,E) ^Wo{E - epf, E-^ep. (23) 

If one takes W^o ^ 2.5°K-^ and^d^'' {pp /dp = Ti-pp/M^, 
My/M = 0.76, from earlier worku3 where a model form of 
W{p, E) was fitted to the measured specific heat in '^He, 
AcBF — 21 is obtained, which is 2 order of magnitude 
larger than both the DMC and the perturbative estimate. 

The effective mass M* is related to the dispersion at 
Pp of the quasi-particle energya 

e{p)^^ + mP,eip)), (24) 

according to h^pp/M* = de{pp)/dp. ThusS'i0 
M*/M = MeMk with the K-mass 

Mk" = 1 + 7?-4-3^S(p, E) |B=e,.p=p, . (25) 
Ti pp op 

and the E-mass 

Me = 1- ^^S(P' E) \E=e^„^ Z-\ (26) 

Thus, the E-mass is nothing but-^^, i.e., the strength 
of the quasi-particle pole at pj-eiJ'EJ. We lack, however. 
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an estimate of the K-mass, to predict M* . Hence, we 
are planning to perform variational and transient esti- 
mate calculations of M* , along the lines of anretpivalcnt 
calculation for the 2-dimensional electron gascScS. 
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FIG. 12. The one-body density matrix p(r) of ^He at 
T = 0°K and pe = 0. 016351^ (full curve). The dashed curve 
gives the product p(i{r)pB{r), with po{r) and ps(r-) respec- 
tively the density matrix for the ideal (uncorrelated) Fermi 
gas and the density matrix of a Bose '^He. 

At the Fermi wave vector pp^ the momentum distribu- 
tion n{p) has a discontinuity and, according to perturba- 
tion theory and to CBF, at least the_additional singular 
behavior of Eq. (^8|). It is knowna that singularities 
dominate the large distance behavior of the Fourier trans- 
form of a generalized function, such as n{p). This imply 
in particular that as r — > oo p{r) ~ ZpQ{r), to leading 
order. On the other hand, using the parameters given in 
Table [II one can show that the first zero of the DMC 



density matrix at the equilibrium density (see, also. Fig. 
^ is at ppf = 5.62, which is not very different from 
ppr — 5.72 — the location of the first zero of po{r). We 
therefore consider the approximate decoupling 



(27) 



with the function pb{t) that must satisfy pb(0) = 1, 
limr^oo Psif) = Z , and we further choose it to be every- 
where non-negative. It is tempting to take psij") as the 
density matrix of a suitable Bose system, with a conden- 
sate fraction Z . 

We have thus simulated a system of '^He with Bose 
statistics, at the equilibrium density of real '^He. We find 
a condensate fraction of 0.208(5), which within error bars 
agrees with our best estimate Z = 0.21(2). In Fig. ^ 
we compare our fit to the density matrix of '^He with 
the prediction of the approximate formula (27), using for 



psir) the fit to the simulated density matrix of a Bose 
■^He. It is apparent that the decoupling of Eq. (^) ap- 
proximately holds. In fact Eq. (^7|) overestimates the 



kinetic energy by 14%, while underestimating the enve- 
lope of the tail of p{r) by about a 20%. Thus statistics 
and correlations in ^He decouple, within a reasonable ac- 
curacy, into those of an ideal (uncorrelated) Fermi gas 
and a Bose ^He, as far as p[r) is concerned. 



V. CONCLUSIONS 

We have presented in this paper results of DMC cal- 
culations of the one-body density matrix and the mo- 
mentum distribution of liquid ^He and liquid '^He. These 
DMC simulations are based on accurate trial wave func- 
tions with fully optimized pair and triplet pseudopo- 
tentials. For ■^He the Fermion sign problem has been 
avoided, by resorting to the fixed-node approximation. 
In particular, we have used backflow nodes, which are 
more accurate than the simpler plane- wave nodes. We 
have recorded our data in a form suitable for future use, 
in terms of analytical fits. 

Our prediction for the condensate fraction of ^He is 
lower than in previous microscopic calculations and we 
have argued that these seems to be consistently related 
to the improved description of "^He afforded by the cal- 
culations presented here. We have also found that our 
results are statistically consistent with the presence of a 
singular term in n{p) as predicted earlier by Gavoret and 
Nozieres, in spite of the absence of explicit long-range 
correlations in our trial functions. 

The discontinuity Z of n[p) at pp in '^He is also pre- 
dicted from our calculations to be sensibly lower than in 
previous variational calculation and in substantial agree- 
ment with GFMC estimates. We have investigate the 
presence in the n(p) of logarithmic singularities, as pre- 
dicted by approximate treatments. We find that our re- 
sults are compatible with the presence of such terms. 
However the strength of such term is in order of mag- 
nitude agreement with perturbation theory, whereas a 
CBF treatment with empirical parameters for the imagi- 
nary part of the optical potential imply a strength which 
differs from our prediction by 2 orders of magnitude. 

We have have also demonstrated that an approximate 
decoupling p(r) ~ po(^)ps ('') holds, with po{r) and 

TABLE IV. Discontinuity of the momentum distribution 
at pf, Z, as function of the density, from various calculations. 
SBF and LBF denote the DMC estimates, obtained using 
the fit of Eq. (^ and results for EMC wave functions with 
short and long-range backflow, respectively. VMC gives the 
variational estimate for the SBF wave function. 



p 


0.01198 


0.01413 


0.01635 


0.01797 


0.01946 






0.348 


0.275 


0.244 


0.221 


VMC 






0.272(2) 






GFMC'' 






£ 0.2 






SBF 


0.46(2) 


0.31(2) 


0.24(1) 


0.14(1) 


0.14(1) 


LBF 






0.21(2) 






''Ref. 
''Ref. 


iB 
l| 
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Pb {r) respectively the density matrix for the ideal Fermi 
gas and the density matrix of a Bose ^He. Thus statistics 
and correlations effects seem to decouple in '^He, as far 
as the density matrix is concerned. The Bose ^He has in 
fact a condensate fraction which agree within error bars 
with the discontinuity Z found in the Fermi ^He. 

We believe that the estimates given in this paper pro- 
vide the most accurate information of this kind on He 
available to date. This is of particular importance for the 
condensate fraction in '^He, as its extraction from deep in- 
elastic neutron scattering seems still not feasible. Using 
our accurate EMC trial functions we have also studied 
partially polarized ^He. We shall report on this study 
elsewhereCJ. 



ACKNOWLEDGMENTS 

Most of the work presented in this paper was done 
when SM was enjoying a postdoctoral fellowship at the 
Laboratorio FORUM of the Istituto Nazionale di Fisica 
della Materia, Pisa. Access to the computing facilities of 
CNUCE are also gratefully acknowledged. 



^ Momentum Distributions, R.N. Silver and P.E.Sokol, eds 
(Plenum, N.Y., 1989). 

See, e.g., A.L. Fetter and J.D. Walecka, Quantum Theory 
of Many-Particle Systems (McGraw-Hill, New York, 1971). 

^ D.M.Ceperley and M.H.Kalos, in Monte Carlo Methods in 
Statistical Physics, K. Binder, ed. (Springer, N.Y., 1979). 
S.Fantoni, Nuovo Cimento 44A, 191 (1978). 

^ M.L.Ristig, in From Particle to Nuclei, Course LXXIX 
Int. School of Phys. Enrico Fermi, A.MoIinari ed., (North- 
Holland, Amsterdam, 1981) 

S.Fantoni and V.R.Pandharipande, Nucl. Phys. A427, 473 
(1984). 

^ E.Manousakis, V.R.Pandharipande and Q.N.Usmani, 

Phys. Rev. B31, 7022 (1985). 
* A. Fabrocini, S. Rosati, Nuovo Cimento Dl, 567, 615 

(1982). 

^ M.F.Flynn, Phys. Rev. B33, 91 (1986). 
A. Fabrocini, V.R.Pandharipande and Q.N.Usmani, Nuovo 
Cimento 14D, 469 (1992). 

E. Feenberg, Theory of Quantum Liquids (Academic, New 
York, 1969). 

Q. N. Usmani, S. Fantoni, V. R. Pandharipande, Phys. Rev. 

B26, 6123 (1982), and references therein. 
" R. A. Smith, Phys. Lett. 63b, 369 (1976); 85B, 183 (1976); 

R. A. Smith, A. Kallio, M. Puoskari, P. Toropainen, Nucl. 

Phys. 328A, 186 (1979). 
" E. Krotscheck, Phys. Rev. B33, 3158 (1986). 

X. Q.G.Wang, S.Fantoni, E.Tosatti and L.Yu, Phys Rev 

B49, 10027, 1994. 

S.Moroni, S.Fantoni, G.Senatore, Europhys. Lett, 30, 93 
(1995). 



" S.Moroni, S.Fantoni, G.Senatore, Phys. Rev., B52, 13547 
(1995). 

^® P.A.Whitlock and R.Panofl, Ca. J. of Phys. 65, 1409 
(1987). 

D.M.Ceperley and E.L.PoUock, Phys. Rev. Lett. 56, 351 
(1986); Can J. of Phys. 65, 1416 (1986); E.L.PoUock and 
D.M.Ceperley, Phys. Rev. B36, 8343 (1986). 
^° G.L.Masserini, L.Reatto, S.A.Vitiello, Phys. Rev. Lett. 69, 
2098 (1992). 

P.J.Reynolds, D.M.Ceperley, B.J. Alder, W.A.Lester, J. 
Chem. Phys. 77, 5593 (1982). 

D. M. Ceperley, in Recent Progress m Many-Body Theo- 
ries, ed. J. Zabolitzky (Springer 1981); D. M. Ceperley, J. 
Stat. Phys. 63, 1237 (1991). 

K.S. Liu, M.H. Kalos, and G.V.Chester, Phys. Rev. AlO, 
303, (1974). 

J. CasuUeras and J. Boronat, Phys. Rev. B 52, 3654, 
(1995). 

S. A. Vitiello, K. E. Schmidt, Phys. Rev. B46, 5442 (1992). 
R. A. Aziz, V. P. S. Nain, J. S. Carley, W. L. Taylor, G. 
T. Conville, J. Chem. Phys. 70, 4330 (1979). 
K. E. Schmidt, M. A. Lee, M. W. Kalos, G. V. Chester, 
Phys. Rev. Lett. 47, 807 (1981). 

E. Manousakis, S. Fantoni, V. R. Pandharipande, Q. N. 
Usmani, Phys. Rev. B28, 3770 (1983). 

2^ Y. Kwon, D. M. Ceperley, R. M. Martin, Phys. Rev. B48, 
12037 (1993). 

^° R.M.Panoff, J.Carlson, Phys. Rev. Lett. 62, 1130 (1989). 
^'^ J.Gavoret and P.Nozieres, Ann. Phys. (NY) 28,349 (1964); 
^2 W.N.Snow,Y.Wang and P.E.Sokol, Europhys. Lett. 19, 403 
(1992) 

P.E.Sokol in Bose-Einstein condensation, A. Griffin, 
D.W.Snoke and S. Stringari, eds. (Cambridge University 
Press, 1995). 

T.R.Sosnick, W.M.Snow, P.E.Sokol and R.N.Silver, Euro- 
phys. Lett. 9, 707 (1989). 

S.Fantoni, B.L.Friman and V.R.Pandharipande, Nucl. 
Phys. A399, 51 (1983). 
^'''R.Sartor and C.Mahaux, Phys. Rev. C21, 1546 (1980). 
P.E.Sokol, K.Skold, D.L.Price, R.Kleb, Phys. Rev. Lett. 54, 
90 (1985). 

S.Fantoni, V.R.Pandharipande and K.E.Schmidt, Phys. 
Rev. Lett. 48, 878 (1982). 

J.P.Jeukenne, A.Lejeunne and C.Mahaux, Phys. Rep. 25, 
83 (1986). 

A.B.Migdal, JETP (Sov. Phys) 5, 333 (1957). 
"1 J.M.Luttinger, Phys. Rev. 119, 1153 (1960). 

Y. Kwon, D.M. Ceperley, and R.M. Martin, Phys. Rev. 
B50, 1684 (1994). 

Y. Kwon, D.M. Ceperley, and R.M. Martin, Phys. Rev. 
B53 7376 (1996). 

M.J. Lighthill, Introduction to Fourier Analysis and Gener- 
alized Functions (Cambridge University Press, Cambridge 
1959). 

'^^ S. Moroni, S.Fantoni, G.Senatore, in preparation. 



10 



